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ABSTRACT 

We report additional detections of the chloronium molecular ion, H 2 CI+, to¬ 
ward four bright submillimeter continuum sources: G29.96, W49N, W51, and 
W3(OH). With the use of the HIFI instrument on the Herschel Space Obser¬ 
vatory, we observed the 2i2 — Iqi transition of ortho-H^^CR at 781.627 GHz in 
absorption toward all four sources. Much of the detected absorption arises in 
diffuse foreground clouds that are unassociated with the background continuum 
sources and in which our best estimates of the A(H 2 GH)/A(H) ratio lie in the 
range (0.9 - 4.8) xl0“®. These chloronium abundances relative to atomic hy¬ 
drogen can exceed the predictions of current astrochemical models by up to a 
factor of 5. Toward W49N, we have also detected the 2i2 — Iqi transition of 
ortho-H 2 ^GH at 780.053 GHz and the In — Oqo transition of para-H^^GH at 
485.418 GHz. These observations imply Hl^GH/H^^GH column density ratios 
that are consistent with the solar system ^^Gl/^^Cl isotopic ratio of 3.1, and 
chloronium ortho-to-para ratios consistent with 3, the ratio of spin statistical 
weights. 

Subject headings: ISM: molecules - Submillimeter: ISM - Molecular processes 


1. Introduction 


The chloronium molecular ion, H2C1’'’, is one of several new interstellar molecules 
discovered using the HIFI instrument (de Graauw et ah 2010) on board the Herschel Space 
Observatory (Pilbratt et ah 2010). Chloronium was hrst detected by Lis et ah (2010) 

- who used the Herschel/B.lFl to perform absorption line spectroscopy of the In — Oqo 
485 GHz transition of para-chloronium toward the bright submillimeter continuum sources 
NGC 66341 and Sgr B2 (S). Toward NGC 66341, the 2i2 — Iqi 782 GHz transition of 
ortho-chloronium was also detected, leading to an inferred ortho-to-para ratio of 3.2, 
close to the ratio of statistical weights for the two spin states. Interstellar H 2 CH was 
subsequently detected by Neufeld et ah 2012 (Paper I) in absorption toward W31C and 
Sgr A (+50kms“^ cloud), and in emission from the Orion Bar photodissociation region 
and Orion South molecular condensation. The chloronium abundances inferred from the 
absorption line detections in Paper I, amounting to as much as 12% of the gas-phase 
chlorine budget, were surprisingly high relative to predictions obtained using the Meudon 
PDR model (Le Petit et ah 2006) with a chemical network for Cl-bearing molecules 
similar to that described by Neufeld & Wolhre (2009, hereafter NW09). These Herschel 
observations were followed by ground-based (IRAM 30-m and CSO) detections of the 
189 GHz lio — Iqi transition of ortho-H 2 Cl”'' toward W31C and W49N (Gerin et ah 2013), 
which - combined with Herschel/Y[\¥\ detections of the In — Oqo transition of para-H 2 GH 

- provided additional estimates of the H 2 GH ortho-to-para ratio. These estimates were 
consistent with the high-temperature LTE value of 3, but had very large error bars. Most 
recently, extragalactic chloronium has been detected in the 2 ; = 0.89 absorber in front of 
the lensed blazar PKS 1830-211 (Muller et al. 2014), again through observations of the 
111 — Oqo transition, redshifted now to ~ 257 GHz. 
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The chemistry of interstellar chloronium has been discussed extensively in NW09 and 
Paper I and only a brief summary is needed here. With an ionization potential (12.97 eV) 
slightly below that of atomic hydrogen, chlorine is singly-ionized in diffuse atomic clouds. In 
regions where the H 2 fraction becomes signihcant, Cl^ reacts exothermically with H 2 to form 
HCl”*", a molecule that has now been detected (de Luca et ah 2012) in the diffuse interstellar 
medium by Herschel/Y{\¥\. A further exothermic reaction with H 2 leads to H 2 CI+, which 
ultimately undergoes dissociative recombination to form HCl or Cl. Thus the overall 
pathway leading to chloronium may be written: Cl (hz/, e) C1^(H2, H) HC1^(H2, H) H2C1’''. 
The chloronium abundance is limited by the dissociative recombination of HCl’*’ and H 2 CI’'’, 
and by the relatively modest flux of UV photons of sufficient energy to initiate the reaction 
sequence by ionizing Cl. Detailed predictions for the depth dependence of the chloronium 
abundance have been presented by NW09, but the predicted column densities were up 
to an order of magnitude lower than the observed values reported in Paper I. Although 
several possible explanations of this discrepancy were investigated in Paper I, none proved 
entirely satisfactory and thus the high chloronium abundances measured toward Sgr A and 
particularly W31C remain puzzling. Like other triatomic hydrides such as H 2 O, H 2 O’*’ and 
H 2 S, chloronium possesses two symmetry states - ortho-H 2 Cl"'’ with total H nuclear spin of 
1, and para-H 2 CD with total spin 0 - each associated with distinct sets of rotational states 
of even (ortho-H 2 Cl’'’) and odd (para-H 2 Cl’*") spatial symmetry. 

Among all the elements lighter than nickel, chlorine possesses the minor isotope (^^Cl) 
that is most abundant relative to its principal isotope (^^Cl), with a solar system ^^Cl/^^Cl 
ratio of only 3.1 (Lodders et ah 2003)111 Thus, the H^^CD isotopologue is often detected 

*^This fact, resulting in an atomic weight of 35.45 for terrestrial chlorine that was unexplained until the 
discovery of isotopes by Aston in 1919, was of historical importance in Turner’s (1833) refutation of Front’s 
Law (which had asserted that the atomic weights of the elements were integral multiples of that of hydrogen). 
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together with the frequencies for both isotopologues falling easily within the 

simultaneous bandwidth of the instrument; indeed, the In — Oqo transition was 

detected in the NGC 63341, W31C, Sgr B2 and PKS 1830-211 observations discussed above. 
In the absence of isotope-selective fractionation in the interstellar gas, the 
ratio serves as a valuable probe of the ^^Cl/^^Cl isotopic ratio outside the solar system; in 
all cases reported to date, including the z = 0.89 absorber in front of blazar PKS 1830-211, 
the ratio is consistent, within the observational uncertainties, with the solar 

system ^^Cl/^^Cl isotopic ratio. 

In this paper, we report new Herschel/Y{W\ observations of interstellar chloronium. 
We have targeted the 2i2 — loi transition of ortho-H 2 CK near 781 GHz, which presents a 
larger optical depth than either of the 485 GHz or 189 GHz transitions observed previously. 
In §2 below, we discuss the observations and data redaction, which have led to nnequivocal 
detections of ortho-H^^GH absorption in foreground diffuse clouds lying along the sight-lines 
to the bright submillimeter continuum sources G29.96, W49N, W51, and W3(OH), together 
with unequivocal detections of ortho-H^^GK and para-H^^GH toward W49N. The resultant 
spectra, and the inferred chloronium column densities and abundances, are presented in §3, 
together with the ortho-to-para ratios and H^^GH/Hg^GH ratios derived for clouds along 
the sight-line to W49N. A discussion of these results follows in §4. 


2. Observations and data reduction 

All the data presented here were obtained as part of the program OTl_dneufeld_l on 
2013 March 26 - 28 in the hnal weeks of the Herschel mission. The 2i2 — loi transitions 
of ortho-H 2 ®GH and H|^GH were observed simultaneously in the lower sideband of Band 
4a of HIFI, and the In — Oqo transitions of para-H^^GH and H|^GH were observed in the 
lower sideband of Band la. These ortho- and para-chloronium transitions were targeted 
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in observations of total dnration 2545 s and 3920 s respectively toward each sonrce, which 
provided spectra with r.m.s. noise levels of ~ 10 niK and ~ 4 mK, respectively. The 
observing time was divided among three separate local oscillator (LO) settings; because 
HIFI employs double sideband receivers, the use of multiple LO settings allows one to 
determine whether any given feature lies in the upper or lower sideband of the receiver. 
Table 1 lists, for each target, the coordinates, distance, Herschel observation identihers 
(OBSID), observed continuum brightness temperatures, and source systemic velocity 
relative to the local standard of rest (LSR). Table 2 lists the rest frequencies for the 2i2 — loi 
and 111 — Ooo transitions of each isotopologue, using data downloaded from the Cologne 
Database for Molecular Spectroscopy (CDMS; Muller et al. 2005), which are based on 
laboratory spectroscopy performed by Araki et al. (2001) and a theoretical estimate of the 
dipole moment by Muller (2008). Because ^^Cl and ^^Cl both possess a nuclear spin of 3/2, 
the 2 i 2 — loi transition is split into eight hyperhne components and the In — Oqo transition 
into three. 

The data were reduced using methods similar to those discussed in Paper I. Because 
of the presence of strong interloper emission lines and/or the small signal-to-noise ratios 
obtained in some of the observed spectra, only six of the 16 possible combinations of target 
line and background source proved to be of value in providing quantitative results: the 
ortho-H 2 ^CD spectra observed toward all four background sources, and the ortho-Hl^CD 
and para-H|^CD spectra obtained toward W49N. In Figures 1 through 3, we show the six 
spectra listed above, with the three LO settings shown separately in red, green and blue. 
Here, the spectra obtained for the two orthogonal polarizations were found to be very 
similar and were averaged to reduce the noise level, yielding signal-to-noise ratios on the 
continuum in the range ~ 100 — 300. Alternate panels show the frequencies computed for 
either the signal (black frames) or image sideband (brown frames) in the rest frame of the 
source. Magenta dashed lines indicate identihed emission lines in the relevant sideband. 


The different spectra obtained at three separate LO settings and in two orthogonal 
polarizations were averaged with a weighting inversely proportional to the square of the 
r.m.s. noise, except for spectral regions where a contaminating emission feature from the 
image sideband caused interference with the H 2 C 1 ^ absorption feature; such regions were 
excluded from the average. In the para-H^^Cl^ spectra observed toward W49N, a ripple 
is clearly apparent in the data and is a known artifact of the instrument (Roelfsema et 
al. 2012). It was removed by using the task fitHiFiFringe within the Herschel Interactive 
Processing Environment version 12.1.0 (HIPE; Ott 2010) to subtract standing waves at 92 
and 98 MHz. Finally, the averaged absorption spectra were converted to a transmission 
fraction under the assumption that the sideband gain is unity; the transmission is then 
given by 2TA/T^(cont) — 1, where Ta is the double-sideband antenna temperature and 
TA(cont) is the continuum temperature. 


3. Results 

In Figure 4, we show (black histograms) the transmission fractions for ortho-H^^CH 
observed toward all four background sources, and for ortho-H^^CH and para-H^^CH 
observed toward W49N. All six panels indicate widespread absorption by foreground 
material unassociated with the background continuum source. We used a multigaussian 
fit to the absorbing column density to deconvolve the hyperfine structure, with the fit 
to the data shown in red and the hyperfine-deconvolved fit in blue. The latter is the 
spectrum that would have resulted in the absence of hyperfine splitting, and shows the same 
frequency-integrated optical depth as the original data. The emission feature appearing at 
LSR velocity ~ 40kms“^ in the G29.96-0.02 ortho-H^^CR spectrum and at LSR velocity 
~ 0kms“^ in the W51 spectrum is a blend of dimethyl ether transitions lying in the same 
sideband as the ortho-H^^CR transition; spectral regions within 10kms“^ of that feature 
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were excluded from the hyperfine-deconvolution procedure. 

In Table 3, we present the column densities derived for ortho-H|^Cl+ within those 
velocity intervals adopted by Indriolo et ah (2015) for each source. These column densities 
were derived for an assumed excitation temperature of 2.73 K, the temperature of the cosmic 
microwave background, an approximation that is expected to be excellent for foreground 
diffuse clouds where the density and radiation held are low, but one that is likely invalid 
for material associated with the background source itself. For an excitation temperature of 
2.73 K, the ratio of the ortho-H 2 Cl+ column density to the total, velocity-integrated optical 
depth for all hyperhne components in the 2i2 — Iqi transition is 4.65 x 10^^ cm“^ perkms“^, 
while the ratio of the para-H 2 Cl+ column density to the total, velocity-integrated optical 
depth for all hyperhne components in the In — Oqo transition is 2.24 x 10^^ cm“^ perkms“^. 

Table 3 also lists the HI column densities given by Indriolo et al. (2015) for each 
velocity interval, based upon recent 21 cm data obtained by Winkel et al. (2015), and the 
corresponding A^(H 2 CH)/A^(H) column density ratios, for an assumed ortho-to-para ratio 
of 3 and an assumed ratio of 3.1. These results are presented graphically in 

Figure 5. Here, dashed horizontal bars indicate velocity intervals close the source systemic 
velocity, where the column density determinations may be unreliable, and vertical error 
bars indicate uncertainties. The latter are highly asymmetric because they include the 
possibility that weak unidentihed emission lines in the relevant spectral window result in an 
underestimate of the inferred H 2 CH optical depth. Based upon an examination of the entire 
spectra shown in Figure 1 - 3, we conservatively adopt a maximum antenna temperature of 
50 mK and linewidth of 5kms“^ for any such emission lines in ortho-chloronium spectra; 
for the para-chloronium spectrum obtained toward W49N, we adopt instead a maximum 
antenna temperature of 25 niK. 

In Table 4, we compare the column densities derived for ortho-H|^CH, ortho-H^^CH, 
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and para-H^^Cl’*' toward W49N. The corresponding ortho-to-para ratios and 
ratios are also tabulated for the velocity intervals adopted by Indriolo et ah (2015). In 
Figure 6, these ratios are plotted in velocity bins of width 5kms“^, with vertical error bars 
again indicating the uncertainties that could result from unidentihed emission lines in the 
spectral window (with a width of 5kms“^ and an assumed maximum antenna temperature 
of 50 mK near 781 GHz and 25 mK near 485 GHz). Values close the systemic velocity of 
the source (8kms“^) are unreliable. 


4. Discussion 

4.1. Chloronium abundances relative to atomic hydrogen 

Our best estimates of the iV(H 2 GH)/iV(H) column density ratios, listed in Table 3 
and plotted in Figure 5, lie in the range (0.87 — 4.84) x 10“®. This range excludes velocity 
intervals close to the systemic velocity of the background source, where the iV(H 2 GH) 
estimates are unreliable. Theoretical predictions, presented previously in Paper I and based 
upon the Meudon PDR code (Le Petit et ah 2006) and the chemical network of Neufeld & 
Wolhre (2009) for the chemistry of halogen-bearing molecules in diffuse clouds can predict 
7V(H2GH)/iV(H) ratios as high as 1 x 10“®, given a gas-phase elemental abundance of 
1 X 10“^ for chlorine (Moomey et al. 2012). Thus, some of the observed iV(H 2 GH)/iV(H) 
ratios exceed the predictions by a factor of at least 5, a discrepancy noted previously 
toward other sources and discussed in Paper I. Moreover, if weak, unidentihed emission 
lines are present within the spectral region where iV(H 2 GH) absorption is observed, the 
actual H 2 GH abundances might signihcantly exceed our best estimates, as indicated by the 
asymmetric error bars in Figure 5. 

None of our best estimates for the V(H2GH)/7V(H) ratio is as large as the value of 
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12 X 10“® reported in Paper I for the sight-line to W31C. However, that large value may 
have been based upon an underestimate of the HI column density toward W31C, which had 
been estimated from 21 cm line spectra presented by Fish et ah (2003). More recent HI 21 
cm observations, carried out by Winkel et ah (2015) using the JVLA, have suggested that 
the 21 cm absorption may in fact be saturated along the sight-line to W31C, providing 
only lower limits on the HI density (Indriolo et ah 2015) and therefore upper limits on the 
7V(H2Cl+)/iV(H) ratio. 

If we now disregard the large iV(H 2 CH)/iV(H) ratio reported for W31C, the maximum 
discrepancy between theory and observation is reduced from the factor ~ 10 discussed 
in Paper I to a factor ~ 5. Moreover, since the publication of Paper I, Novotny et ah 
(2013) have reported a laboratory measurement of the rate coefficient for dissociative 
recombination (DR) of HCH, a key intermediary in the primary production pathway for 
chloronium: Cl (hz/, e) CH(H 2 , H) HCH(H 2 , H) H 2 CD. The measured value is smaller than 
that adopted by NW09 by a factor of 1.5 at 10 K and 3 at 300 K. Since the DR of HCD 
limits the availability of HCD to form H 2 CD, a downward revision in the adopted rate 
coefficient would increase the predicted H 2 CD abundances; however, the dependence is 
a relatively weak one, in that a factor 10 reduction in that rate coefficient was found to 
increase the H 2 CD abundance by only a factor of 1.3 (Paper I). Definitive measurements 
of the DR rates for HCD under astrophysical conditions - and also for H 2 CD, which is 
primarily destroyed by DR - will have to await future experiments at the Cryogenic Storage 
Ring (CSR) in Heidelberg. This facility will permit the hrst storage ring measurements 
of rotationally-cold molecular ions, in contrast to previous measurements in which the 
recombining ions had a rotational temperature ~ 300 K. The case of rotationally-cold ions 
is most appropriate to hydride molecular ions in low-density diffuse clouds, where most ions 
are in the ground rotational states; depending upon the exact details specific to each DR 
process, future CSR measurements could yield rate coefficients that are either smaller or 
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larger than those obtained with rotationally-warm ions. 

The astrocheniical models presented in Paper I predict that the maximum 
A^(H2C1+)/A^(H) ratios are attained in diffuse clouds of total visual extinction 0.3 mag) 
in which the gas is primarily atomic, a prediction that is supported by the observed 
distribution of H2C1'^ in velocity space. In particular, the H2C1^ absorption line profiles 
plotted in Figure 1 are more similar to those of OH”*" and H20^ (Indriolo et al. 2015) 

- two species believed to originate in gas of small molecular fraction - than they are to 
those of H 2 O, HF (Sonnentrucker et al. 2015) or CS (Neufeld et al. 2015). The latter three 
molecules, believed to be most abundant in gas with a high molecular fraction, tend to show 
narrow absorption features that are entirely absent in the OH+, H20^ and H2C1”'' spectra. 


4.2. Chloronium isotopologue ratio 

Because isotope-selective fractionation is expected to be entirely negligible, the 
A^(H 2 ^Cl'^)/iV(H 2 ^Cl''") ratio provides a valuable measure of the ^^Cl/^^Cl isotopic ratio 
within the interstellar gas. This probe has been exploited previously by Lis et al. (2010) 
and in Paper I, with observations of para-chloronium absorption implying best-fit ^^Cl/^^Cl 
ratios in the range ~ 2 — 4 within foreground clouds along the sight-lines to NGC 63341, 
W31C and Sgr A (-1-50kms“^ cloud); these values bracket the solar system ^^Cl/^'^Cl 
ratio of 3.1 (Fodders 2003). Similar results have been obtained from observations of HCl 
emission toward W3 A (Cernicharo et al. 2010) and several other Galactic sources (Peng et 
al. 2010). In the one extragalactic determination of a A^(H 2 ^GG)/A^(H 2 ^GG) ratio obtained 
to date, Muller et al. (2014) reported a value of 3.11 q; 2 for the 2 ; = 0.89 absorber in front 
of PKS 1830-211, based on absorption-line observations of para-chloronium. In the present 
study, our observations of the stronger ortho-chloronium transitions towards W49N permit 
more accurate determinations of the ^^Gl/^^Gl ratio than were possible in Paper I. As 
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shown in Figure 6 (middle panel), the iV(H 2 ®Cl’^)/iV(H 2 ^Cl’'') ratios determined in 5kms“^ 
velocity bins are all consistent with the solar system ^®C1/^^C1 isotopic ratio to within the 
uncertainties. The average for the entire sight-line is 3.50lgg2 (Table 4). 

The solar system ^®C1/^^C1 ratio is not entirely understood. Citing Galactic chemical 
evolution models of Kobayashi, Karakas & Umeda (2011; hereafter KKUll), Muller et al. 
(2014) noted that the isotopic ratios for most elements in the 2 ; = 0.89 absorber in front 
of PKS 1830-211 were htted well by a KKUll model for stars in the solar neighborhood 
with [Fe/H] = -2.6 (i.e. stars with iron abundances, relative to H, a factor 10“^'® times 
the solar value). For that model, in which Type II supernovae were the source of heavy 
elements, the predicted ^^Cl/^^Cl isotopic ratio was 3.0, very close to the solar system value 
of 3.1 and entirely consistent with the value of 3.11 q 2 for the z = 0.89 absorber in front of 
PKS 1830-211. However, the KKUll model with [Fe/H] = 0.0 (i.e. with solar abundances) 
predicts a ^^Cl/^^Cl isotopic ratio of only 1.9, a factor ~ 1.6 below the solar system value, 
and a value that is signihcantly discrepant with the results we have obtained here for the 
Galactic ISM. This modest discrepancy notwithstanding, the KKUll model does predict 
that the ^®C1/^^C1 will show a relatively weak dependence on [Fe/H], in sharp contrast to 
most of the other elemental isotopic ratios considered. Thus, although the solar ^^Cl/^^Cl 
ratio is evidently not predicted with perfect accuracy, the model does correctly suggest that 
35 qi/37qi evolve strongly. 


4.3. Chloronium ortho-to-para ratio 

The ortho-to-para ratios (OPRs) plotted in Figure 6 (bottom panel) are all consistent 
with the ratio of statistical weights, 3, to within the estimated uncertainties. This is the 
behavior expected in LTE at temperatures much larger than AE^p/k = 20.2 K, where AEop 
is the energy difference between the lowest rotational states of ortho- and para-chloronium 
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(loi and Ooo respectively.) Observations of some other polyatomic hydrides in diffuse 
molecular clouds reveal OPRs consistent with the ratio of statistical weights - examples 
include H 2 O (Flagey et ah 2013; Lis et ah 2013) and H 20 ’'' (Schilke et ah 2013; Gerin et 
ah 2013) - but this behavior is not universal. In particular, Hj]' (Crabtree et ah 2011) and 
NH 3 (e.g. Persson et ah 2012) both show OPRs with signihcant departures from the values 
expected in LTE, a behavior that is believed to reflect the OPR of H 2 molecules involved 
in the formation of these hydrides through nuclear-spin-conserving gas-phase reactions (e.g. 
Oka 2004, Faure et ah 2013; Le Gal et al. 2014). 


When chloronium is formed by the reaction HCR(H 2 , H)H 2 CR - or indeed when 
any triatomic hydride molecular ion (such as H 20 ^) is formed by an analogous hydrogen 
abstraction reaction - the resultant ions are formed with an initial ortho-to-para ratio, 
OPRo, that is determined by the nature of the intermediate complex; if the latter is 
long-lived, then the reaction occurs by the scrambling mechanism (e.g. Oka 2004, Herbst 
2015), in which all routes that lead to products are allowed. In this case, one-half of 
all reactions of HCR with para-H 2 lead to ortho-H 2 Cl’'' (with the remainder forming 
para-H 2 CR), and hve-sixths of all reactions of HCR with ortho-H 2 lead to ortho-H 2 CR. 
Thus, H 2 CR is formed with an initial ortho-to-para ratio, OPRq, that is related to the 


OPR of Ho: 


OPRo(H2C1+) = 


5 0PR(H2) + 3 


0PR(H2) + 3 ■ 

However, if the formation mechanism does not involve a long-lived intermediate complex, 
but occurs instead by a hopping mechanism in which a hydrogen atom breaks off from 
H 2 and attaches to HGH in a long-range interaction, then OPRo(H 2 GR) is the ratio of 
statistical weights (i.e. 3), regardless of OPR(H 2 ) . Detailed quantal calculations will be 
needed to determine which case applies to H 2 GR. 


( 1 ) 


As discussed by Gerin et al. (2013) and Herbst (2015), the relationship between 
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OPRo(H2 C1''') and the observed OPR for H 2 CR depends upon the relative rates of two 
processes: (1) destruction, which occurs primarily via dissociative recombination at a rate 
that we assume to be identical for the two spin symmetries; and (2) thermalization via 
reactive collisions: 

ortho - H 2 CI+ + H ^ para - H 2 CI+ + H. (2) 

The latter process tends to drive OPR(H 2 CR) to its LTE value, OPRlte(H 2 CR). The solid 
black curve in Figure 7 shows the dependence of OPRlte(H 2 CR) upon the gas temperature, 
while the red curve shows the analogous quantity for molecular hydrogen, OPRlte(H 2 ). 
Dashed red and black curves show the functions 9 exp(—20.2K/T) and 9 exp(—170K/T) 
applying in the low temperature limit where only the lowest two rotational states are 
populated. 


Given an H 2 OPR that is in LTE - the behavior typically predicted in theoretical 
models for diffuse molecular clouds - and if the formation of H 2 CD occurs via the scrambling 
mechanism, then the blue curve shows OPRo(H 2 CR), as computed from equation (1). 
The observed OPR(H 2 CR) is therefore expected to lie between the blue and the black 
curves, depending upon the relative rates of destruction and thermalization. Typically, 
astrochemical models (e.g. NW09) predict that chloronium is most abundant in gas that is 
primarily atomic and has a temperature ~ 100 K. 


In steady-state, and with the assumptions given above, we find that 


OPR(HjCl+) = + fcd.».OPRo(H;CC) 

^op^(ff) T ^dr^e 


( 3 ) 


where kgp is the rate coefficient for the forward reaction in (2), kpo is the rate coefficient for 
the reverse reaction, and k^r is the rate coefficient for dissociative recombination of H 2 CD. 
Now using the principle of detailed balance to infer that kpo = /copOPRlte(H2 C1^)5 we may 
obtain the chloronium OPR as a weighted mean of OPRlte(H2 C1^) (black curve in Figure 
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7) and OPRo(H 2 Cl+) (blue curve in Figure 7) 

0 PR(H2C1+) = a; 0 PRlte(H2 C1+) + (1 - x)OPRo(H 2 Cl+), (4) 

where x = ko^/[ko^n{YV) + /cdr’^e] lies between zero and unity. For 100 K gas that is 
primarily atomic with an electron abundance equal to the gas-phase elemental abundance 
of carbon (1.6 x 10“^), we hnd that fcdrne/n(H) ~ 5 x 10“^^ cm^ s“^. Thus, if thermalization 
proceeds at the Langevin rate with kop ~ 10“®, then x in equation (4) is roughly 0.95 
and the expected chloronium OPR is close to OPRlte; in the opposite limit where 
kop -C 5 X 10“^^cm^s“^, we hnd that a; <C 1 and the OPR is close to OPRq. Once again, 
future studies will be needed to estimate the value of kop- 

Given our assumptions about the possible errors, all of the observed chloronium OPR 
values plotted in Figure 5 (lower panel) are consistent with 3, and most are consistent 
with 2.5. In light of the foregoing discussion, and given a strong prior from astrochemical 
modeling that chloronium in diffuse clouds is primarily present in gas with a temperature 
of at least 30 K (for which OPRlte(H2 C1^) is very close to 3), an OPR between 2.5 and 3 
can be explained in any of the following scenarios: (1) chloronium is formed by the hopping 
mechanism; (2) the rate coefficient kop ^ 1.5 x 10“^°cm^s“^ (such that x ^ 0.75); or (3) the 
observed chloronium lies mainly in gas with T ^ 110 K (such that OPRlte(H 2 ) > 1.8 and 
OPRlte(H2 C1^) ^ 2.5). Further work will be needed to determine which of these scenarios 
applies (or apply). 

HIFI has been designed and built by a consortium of institutes and university 
departments from across Europe, Canada and the United States under the leadership of 
SRON Netherlands Institute for Space Research, Groningen, The Netherlands and with 
major contributions from Germany, France and the US. Consortium members are: Canada: 
CSA, U. Waterloo; France: CESR, LAB, LERMA, IRAM; Germany: KOSMA, MPIfR, 
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MPS; Ireland, NUI Maynooth; Italy: ASI, IFSI-INAF, Osservatorio Astrofisico di Arcetri- 
INAF; Netherlands: SRON, TUD; Poland: CAMK, CBK; Spain: Observatorio Astrononiico 
Nacional (IGN), Centro de Astrobiologia (CSIC-INTA). Sweden: Chalmers University of 
Technology - MC2, RSS & CARD; Onsala Space Observatory; Swedish National Space 
Board, Stockholm University - Stockholm Observatory; Switzerland: ETH Zurich, FHNW; 
USA: Caltech, JPL, NHSC. 
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and AYA2012-32032. 
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Table 1: List of observed sources 


Source 

R.A. 

Dec 

Distance® 

71 ^ 

^sys 

OBSID 

Continuum rA(DSB) 


J2000 

J2000 

kpc 

km s“^ 

-1342268000 

(485 GHz) 

(781 GHz) 

G29.96-0.02 

18:46:03.80 

-02:39:22.0 

5.26 

98 

615-617 


1.17 K 






459-461 

0.25 K 


W49N 

19:10:13.20 

+09:06:12.0 

11.11 

11 

618-620 


4.80 K 






456-458 

1.24 K 


W51 

19:23:43.90 

+14:30:30.5 

5.1 

58 

621-623 


5.40 K 






453-455 

1.44 K 


W3(OH) 

02:27:03.82 

+61:52:24.6 

2.04 

-44 

625-627 


1.88 K 






470-472 

0.49 K 



References for distance and Vsys estimates in the order listed: Zhang et al. (2014); Zhang et al. (2013); 
Sato et al. (2010); Hachisuka et al. (2006). All distance estimates are trigonometric estimates based on 
observations of water masers. 






Table 2: Chloronium line frequencies and spontaneous radiative rates 


Species 

Line 



Frequency^ 

(GHz) 

Aui 

(s-^) 

o - Hf C1+ 

2 i2 ~ lot 

F = 5/2 

-5/2 

781.60886 

1.79 X 10-3 



F = 3/2 

-1/2 

781.61075 

2.48 X 10-3 



F = 3/2 

-5/2 

781.62125 

2.98 X 10-^ 



F = 5/2 

-3/2 

781.62240 

4.17 X 10-3 



F = 7/2 

-5/2 

781.62653 

5.96 X 10-3 



F = 1/2 

-1/2 

781.62820 

4.96 X 10-3 



F = 3/2 

-3/2 

781.63479 

3.18 X 10-3 



F = 1/2 

-3/2 

781.65224 

9.93 X 10-^ 

o 

1 

toco 

Q 

2 i 2 — lot 

F = 5/2 

-5/2 

780.03907 

1.78 X 10-3 



F = 3/2 

-1/2 

780.04058 

2.47 X 10-3 



F = 3/2 

-5/2 

780.04886 

2.96 X 10-^ 



F = 5/2 

-3/2 

780.04977 

4.14 X 10-3 



F = 7/2 

-5/2 

780.05305 

5.92 X 10-3 



F = 1/2 

-1/2 

780.05437 

4.93 X 10-3 



F = 3/2 

-3/2 

780.05956 

3.16 X 10-3 



F = 1/2 

-3/2 

780.07335 

9.87 X 10-^ 

p - Hf C1+ 

1 

o 

o 

O 

F = 3/2 

-3/2 

485.41343 

1.59 X 10-3 



F = 5/2 

-3/2 

485.41767 

1.59 X 10-3 



F = 1/2 

-3/2 

485.42080 

1.59 X 10-3 


“Based on laboratory spectroscopy of Araki (2001). Downloaded from the CDMS database on 2015 Mar 25. 
^For an assumed dipole moment /i = 1.89 Debye (Muller 2008). The spontaneous radiative rates are propor¬ 
tional to . 
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Table 3: Chloronium column densities and abundances 


Source 

^^LSR 
(km s“^) 

N(o - h|'*C 1+) 
(10^^ cm“^) 

iV(H2Cl+)“ 
(IqI^ cm“^) 

AAr(H2Cl+)'’ 

(10^^ cm“^) 

JV(H)'= 
(10^^ cm“^) 

JV(H2C1+)/N(H) 

X 10® 

G29.96-0.02 

[-10,6] 

1.29 

2.28 

2.62 

1.40 

1.63 


[6, 17] 

2.68 

4.73 

2.62 

3.45 

1.37 


[17, 28] 

0.99 

1.74 

2.62 

1.42 

1.22 


[56, 65] 

3.85 

6.79 

2.62 

1.62 

4.19 


[65, 73] 

4.33 

7.63 

2.62 

2.32 

3.29 


[73, 79] 

1.18 

2.07 

2.62 

1.69 

1.23 


[79, 88] 

1.49 

2.63 

2.62 

2.58 

1.02 


[88, 95]*^ 

2.00 

3.52 

2.62 

1.84 

1.91*^ 


[95, 113]'^ 

11.45 

20.19 

2.62 

6.02 

3.35'^ 

W49N 

[1, lO]"^ 

7.16 

12.63 

0.64 

4.44 

2.84'^ 


[10, IT]*^ 

12.79 

22.55 

0.64 

3.58 

6.30'^ 


[17, 25] 

4.55 

8.03 

0.64 

2.34 

3.43 


[25,43] 

15.88 

28.01 

0.64 

5.58 

5.02 


[43, 51] 

3.35 

5.91 

0.64 

1.67 

3.54 


[51, 66] 

12.07 

21.29 

0.64 

6.15 

3.46 


[66, 80] 

3.68 

6.48 

0.64 

2.89 

2.24 

W51 

[11,16] 

1.49 

2.63 

0.57 

0.84 

3.13 


[16, 21] 

0.45 

0.79 

0.57 

0.88 

0.90 


[21, 33] 

1.65 

2.91 

0.57 

1.42 

2.05 


[33,42] 

1.02 

1.79 

0.57 

0.84 

2.13 


[42, 55]'^ 

7.27 

12.82 

0.57 

4.43 

2.89'^ 


[55, 62]'^ 

3.29 

5.80 

0.57 

3.25 

1.79'^ 


[62, 75] 

11.85 

20.90 

0.57 

> 1.54 

< 13.57 

W30H 

[-51, -SO]*^ 

8.06 

14.22 

1.63 

1.30 

10.94“^ 


[-25, -8] 

2.63 

4.63 

1.63 

1.39 

3.33 


1 

00 

2.41 

4.25 

1.63 

0.95 

4.47 


“Total chloronium column density (both symmetry states and both isotopologues), for an assumed 
A^(H2^C1+)/N(H2^C1+) ratio of 3.1 and an assumed OPR of 3 

^Amount by which A^(H2C1“*') is underestimated if a weak emission line of integrated intensity 250 mK kms“^ 
is present in the velocity interval 

“From Indriolo et al. (2015), based on HI 21 cm observations by B. Winkel et al. (2015, in preparation) 
‘^Velocity intervals close to the source systemic velocity, for which the 1V(H2C1''') estimates are unreliable 
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Table 4: 

Chloronium column densities, isotopic ratios and ortho-to-para ratios toward W49N 

VLSR 

N{o - Hf C1+) 

N{o- Hf C1+) 

lV(p - C1+) 

iV(HfCl+)/N(HfCl+) 

OPR 

(km 

) ( 10 ^^cm“^) 

( 10 ^^ cm~^) 

( 10 ^^ cm“^) 

Note (a) 

Note (a) 

[ 1 , 10 ]' 

7.16 

2.36 

0.64 

3 04+0-21 

11 16^'^'''’^ 
11-10-4.62 

[10,17]'’ 

12.79 

3.62 

2.69 

q cq+O.lS 
o.OO_o_42 

4-76;°:^« 

[17,25] 

4.55 

0.67 

1.15 

6.8i;0;P 

3-98;?+ 

[25,43] 

15.88 

4.06 

4.55 

4 qi+012 
'^•^-‘■-0.42 

4 40 + 0.11 

[43,51] 

3.35 

1.16 

1.20 

2 SQ+0-42 

2.80;[{;^? 

[51,66] 

12.07 

4.06 

3.79 

2 97+0.12 
^•^'-0.32 

4 TO+O.IS 

[66,80] 

3.68 

1.36 

1.23 

o 7(4+0.36 

2 qq+0.39 

[17,80] 

39.54 

11.31 

11.91 


q qr) + 0.20 

'^•'^^-0.53 


“Error bounds indicate the errors that could result from the presence of a weak interloper emission line within 
the velocity interval of integrated intensity 250 mK kms“^ (125 mK kms“^ at the p — H^^Cl’*" frequency). 
^Velocity intervals close to the source systemic velocity, for which the 7V(H2C1’'') estimates are unreliable 
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G29.96-0.02 781.627 GHz: Signal SB 




W51 781.627 GHz: Signal SB 




Fig. 1.— Spectra for 2i2 — loi transition of ortho-H^^Cl’*' obtained at three LO settings are 
shown separately in red, green and blue, with the frequency scale computed for either the 
signal or image sideband (alternate panels) in the rest frame of the source. Magenta dashed 
lines indicate identified spectral lines in the relevant sideband. Top two panels: G29.96 - 
0.02, for an assumed source velocity of +98 kms“^ relative to the LSR. Bottom two panels: 
W51, for an assumed source velocity of +55 km s“^ relative to the LSR. 































26 


W3(0H) 781.627 GHz: Signal SB 



W3(0H) 781.627 GHz: Imoge SB 



W49N 781.627 GHz; Imoge SB 



Fig. 2.— Same as Figure 1, except for W3(OH) (top two panels, for an assumed source 
velocity of -45 km s“^ relative to the LSR) and W49N (bottom two panels, for an assumed 
source velocity of +8 km s“^ relative to the LSR). 
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W49N 780.053 GHz: Signal SB 



W49N 780.053 GHz: Image SB 





Fig. 3.— Same as Figure 1, except for the ortho-H^^Cl’*’ 2i2 — Iqi (top two panels) and 
para-Hl^Cl'*' In — Oqo transitions (bottom two panels) observed toward W49N. 
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G29,96-0.02 o-H^Cr 781.627 GHz 


W49N o-H, Gr 781.627 GHz 




0.5 F.. I I 

-20 0 20 40 60 80 

LSR velocity (km s”^) 


W51 o-H,^^cr 781.627 GHz 


W49N o-H/xr 780.053 GHz 



W3(OH) o-H 2 ^^cr 781.627 GHz 



-20 0 20 40 60 80 

LSR velocity (km s”'') 

W49N p-H 2 ^^CI^ 485.418 GHz 



1 .1 

1.0 


: iW V ^ j 

0.9 

- \aa/ ^ 

l 1 

g 

1/ V 



w 
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1 1 “ 

cq 
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00 


llfl 

CO 
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1 

0.7 


1 

0.6 
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-60 -40 -20 0 

LSR velocity (km s"') 


0 20 40 60 

LSR velocity (km s”') 


Fig. 4.— Transmission fractions for the ortho-H^^Cl”*' 2i2 — loi transition observed toward 
G29.96 - 0.02, W51, W3(OH) and W49N; and for the ortho-H^^CG 2i2 — loi and para- 
111 — Ooo transitions observed toward W49N. Black histogram: observed data. Red 
curve: £t to data. Blue curve: hyperhne-deconvolved fit. The emission lines in the G29.96 
~ 0.02 and W51 panels are due to dimethyl ether. 
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G29.96-0.02 
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Fig. 5.— Inferred iV(H 2 Cl^)/iV(H) ratios for the velocity intervals adopted by Indriolo 
et al. (2015). Dashed horizontal bars indicate velocity intervals close the source systemic 
velocity, where the column density determinations are unreliable. Vertical error bars indicate 
uncertainties, and include the possibility that weak emission lines (with a width of 5 km s“^ 
an assumed maximum antenna temperature of 50 mK) result in an underestimate of the 
inferred H 2 CD optical depth. 
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10 



o f I 

0 20 40 60 80 

LSR velocity (km s“^) 



0 20 40 60 80 

LSR velocity (km s"^) 


Fig. 6.— Top panel; derived column densities for ortho-H^^Cl’*’ (black), ortho-Hl^Cl"^ (blue), 
and para-H^^Cl’*' (red) along the sight-line to W49N. Values close the systematic velocity of 
the source (8 km s“^) are unreliable. Middle panel: inferred ortho-H^^Cl’*’ / ortho-H|^Cl’'' 
ratio in velocity bins of width 5 km s“^. Bottom panel: inferred ortho-to-para ratio for 
fl|^Cl+ in velocity bins of width 5 km s“^. Vertical error bars in the bottom two panels 
indicate uncertainties, and include the possibility of weak emission lines (with an assumed 
width of 5 km s“^ and maximum antenna temperature of 50 mK near 781 GHz or 25 mK 
near 485 GHz) resulting in an underestimate of the inferred optical depths. 
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0 100 200 300 

Gas temperature (K) 


Fig. 7 .— Ortho-to-para ratios (OPR) as a function of gas temperature. Solid black curve; 
H 2 CR OPR in LTE. Dashed black curve: 9 exp(—20.2K/T). Solid red curve: H 2 OPR in 
LTE. Dashed red curve: 9 exp(—170K/T). Blue curve: initial OPR of H 2 CD formed by 
reaction of HOD with H 2 , for the scrambling mechanism (see text) and with the H 2 OPR in 
LTE. 









